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Abstract. 

Ultra-Relativistic Heavy Ion Collisions at an energy y^sjvjv = 65 GeV 
are studied in a three-dimensional Fluid Dynamical model. The results of a 
hydrodynamical evolution using the PIC-method are shown. The importance and 
diagnostic value of a proper reaction plane determination is emphasized, and the 
time development of collective obscrvablcs is presented. 
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Introduction: Fluid Dynamical (FD) models arc widely used to describe ultra- 
relativistic heavy ion collisions. Their advantage is that one can vary flexibly the 
Equation of State (EoS) of the matter and test its consequences on the reaction 
dynamics and the outcome. In energetic collisions of large heavy ions, especially 
if QGP is formed in the collision, one-fluid dynamics is a valid and good description 
for the intermediate stages of the reaction. Here, interactions are strong and frequent, 
so that other models have limited validity. On the other hand, the initial and final, 
Freeze-Out (FO), stages of the reaction are outside the domain of applicability of the 
fluid dynamical model. 

The equations of perfect FD are just the local energy, momentum and baryon 
number conservation laws, where we assume local equilibrium, so that the energy- 
momentum tensor is given by the EoS. Thus, the most important and almost only 
input is the EoS. 



| Speaker 



Fluid Dynamics as Diagnostic Tool 



2 



Of course the initial and final (F0-) conditions are also vital, and the deviations 
from thermal equilibrium (viscosity and dissipation) are also important. Here, though, 
we will concentrate on the impact of the initial condition. 

As pointed out above, the equations of perfect hydrodynamics are the 
conservation laws. Still, the approach additionally assumes a local thermal 
equilibrium. If this is present, intensive macroscopic space-dependent quantities such 
as temperature, T, chemical potential, fiB, and pressure, P, can be defined. An 
Equation of State has to be introduced, which gives a relationship between energy- 
and baryon number density and the pressure: P = P (e, n). This is necessary to have 
a complete, solvable set of equations. After introducing the shear viscosity, t], and the 
heat conductivity, n, a state near equilibrium can also be described. 

The PIC Model: The Particle in Cell (PIC) method was developed by Amsden and 
Harlow in the early 1960s [1] and upgraded to ultra-relativistic energies by Nix and 
Strottman in the 80s and 90s [2]. The method takes advantage of a combination of 
Eulerian and Lagrangian solution methods; the pressure is calculated on an Eulerian 
grid, while all current transfers are calculated in a larger number of Lagrangian cells, 
the so-called "marker particles" . 

The method has numerous advantages (e.g. keeping conserved quantities 
constant), and a typical main disadvantage, the development of "ringing instabilities" 
in rapid expansion or explosion. This has recently been eliminated by assigning the 
marker particles random positions to avoid structures, along which an instability would 
start to develop. With this upgrade the stability of expansion is increased, so that a 
spherical object remains exactly spherical (within the numerical resolution) even after 
a volume increase by a factor of more than 100. 

After upgrades, the code can follow the fluid dynamical expansion even well 
beyond the physical applicability of FD. This extended FD history can then be used to 
analyse the effects of choices on the FO hyper-surface, or on the FO layer (a 4-volumc 
domain) . 

A code with finite cell-size exhibits necessarily a numerical viscosity and coarse 
graining. The kinetic energy lost in this averaging is thermalized just as in nature. 
No additional constraint is enforced to conserve the entropy. 

The numerical viscosities are in the range of the estimated physical viscosity [3], 
if the reaction is not too peripheral, in which case the applicability of FD is not the 
best anyway. The larger than realistic viscosities in the transverse direction dissipate 
more transverse flow to heat, so we get higher T and higher thermal smearing. 

The initial State: The initial shape is calculated in an effective ID Yang-Mills 
coherent field model, where the transverse plane is split up to "streaks" like in the 
"Firestreak" -model of the 1970s [HE]. Projectile and target matter interpenetrate 
each other initially in the collision, exchange color charge and form "strings" which 
expand until all initial kinetic energy is converted to string tension. No transverse 
momentum exchange is taken into account among the streaks in this initial stage, but 
a longitudinal expansion at the edges of the strings was included in the expansion to 
make the surface at the initial moment (when the FD starts) smoother. 

The streaks expand until the initial kinetic energy is converted fully into string 
tension. We use an effective string tension of 6 — 8GeV/fm, which was deduced 
from earlier string model fits to experimental CERN - SPS results. This results in 
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Figure 1. A schematic view on the initial state of a heavy ion collision as used 
in PIC. The labeled lengths indicate the dependence on the particular quantity, 
not the quantity itself. 



moderately long streaks of 4 — 8 fm, so that the thickness of the initial state is rather 
a flat than an extremely long object parallel to the beam direction. 

While in the middle of the transverse plane in the CM frame the matter is at 
rest, the edges on the projectile and target side have projectile- and target-directed 
excess momenta, which lead to a skewed object at the end of the initial state. The 
longer the formation of initial state lasts, the more skewed this initial state becomes 



Non-existing Symmetry and Event-plane reconstruction: Although both nuclei 
involved in the collision of two gold beams are of equal size and mass, there is 
no forward-backward (F/B)-symmetry. In each event a clear distinction between 
projectile and target hemisphere must be made. It is trivial to see that a single- 
event particle distribution following from an initial state of finite impact parameter 
as described above, is not F/B-symmetric but deviates from this by a certain angle 
Q in the x/z-plane (see figure [T]). If the event-plane reconstruction is not done for its 
F/B-orientation, these effects are lost. 

The aforementioned effect, the third flow component, has been predicted as a 
signature of the QGP [HIE]. It should be visible in the special behaviour of the directed 
flow coefficient v\ (y) , which is supposed to be closer to zero at central rapidities than 
the usual linear increase of v\ vs. rapidity, y. This effect is measurable at around 
midrapidity [5] 

Especially correlation measurements, but also flow determinations, must be aware 
of this fact. Two-particle correlations showing basically two peaks at A?/ = and 
Acj) = 0, 2tt do not take the real geometry into account. Taking a sample over events 
with randomly identified F/B event planes obscures the advent of interesting physics. 

Results: The distorted initial state has the biggest pressure gradient along the lines 
indicated in figure [T] This leads to a distribution of matter and energy density in an 
ellipsoid shape in very late stages of the collision, which is enhanced in the direction of 
the original pressure gradient. Figure [2] shows the calculated energy-density in the x/z- 
plane after 8 (a) and 71 (b) cycles, respectively (corresponding to times t\ — 0.3 fm/c 
and to, — 2.7 fm/c after the beginning of the hydro-evolution). The flow patterns that 
follow from the initial state can be seen in figure O The presence of the third flow 
component is well visible in the former. It can also be seen that both directed and 
elliptic flow coefficients approximately double their magnitude between the two times. 



(figure [J) . 
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Figure 2. (Color online) The energy-density in the y = — plane at two 

different times (after 0.3 fm/c (a) and 2.7 fm/c (b)) in Au+Au-collisions at 

v / S]vjv = 65 GeV at impact parameter b = 0.7 (Rl + R2). Note different scales 
between the plots. 




Figure 3. (Color online) The flow coefficients v\ and V2 as a function of 
pseudorapidity r) at two different times (after 0.3 fm/c (a) and 2.7 fm/c (b)) 
in Au+Au-collisions at ^/sjvjv = 65 GeV at impact parameter b = 0.7 (Ri + Raj- 



All flow plots are calculated assuming massless quarks in a QGP-Phase. In 
evaluating the flow during the hydro evolution we take advantage of the observed 
iVq-scaling. It should further be noted that all flow plots are calculated at constant 
time hypersurfaces. This is of course not including freeze-out and hadronization, but 
it allows us to study effects of different initial conditions during the hydro evolution. 
This leads to a larger v\ than in experimental data [8], because, we neither performed 
a FO calculation through a realistic FO hypersurface, nor an impact parameter 
averaging. Both these processes make the distributions smoother and decrease v\ and 
v 2 ■ The "dip" in w 2 at midrapidity is a consequence of too small numerical viscosity 
in longitudinal direction. 
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